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In  the  design  and  study  of  solar  energy,  information  on  solar  radiation  and  its  components  at  a  given 
location  is  very  essential.  Solar  radiation  data  are  required  by  solar  engineers,  architects,  agriculturists 
and  hydrologists  for  many  applications  such  as  solar  heating,  cooking,  drying  and  interior  illumination  of 
buildings.  For  this  purpose,  in  the  past,  several  empirical  correlations  have  been  developed  in  order  to 
estimate  the  solar  radiation  around  the  world.  The  main  objective  of  this  study  is  to  review  the  global 
solar  radiation  models  available  in  the  literature.  There  are  several  formulate  which  relate  global 
radiation  to  other  climatic  parameters  such  as  sunshine  hours,  relative  humidity  and  maximum 
temperature.  The  most  commonly  used  parameter  for  estimating  global  solar  radiation  is  sunshine 
duration.  Sunshine  duration  can  be  easily  and  reliably  measured  and  data  are  widely  available. 
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1.  Introduction 

Information  of  local  solar  radiation  is  essential  for  many 
applications,  including  architectural  design,  solar  energy  systems 
and  particularly  for  design.  Unfortunately,  for  many  developing 
countries,  solar  radiation  measurements  are  not  easily  available 
due  to  the  cost  and  maintenance  and  calibration  requirements  of 
the  measuring  equipment.  Therefore,  it  is  important  to  elaborate 
methods  to  estimate  the  solar  radiation  based  on  readily  available 
meteorological  data. 

Several  empirical  models  have  been  developed  to  calculate 
global  solar  radiation  using  various  climatic  parameters.  These 
parameters  include  extraterrestrial  radiation,  sunshine  hours, 
mean  temperature,  maximum  temperature,  soil  temperature, 
relative  humidity,  number  of  rainy  days,  altitude,  latitude,  total 
precipitation,  cloudiness  and  evaporation.  The  most  commonly 
used  parameter  for  estimating  global  solar  radiation  is  sunshine 


*  Fax:  +90  442  2360957. 

E-mail  address:  abakird@atauni.edu.tr. 

1364-0321/$  -  see  front  matter  ©  2009  Elsevier  Ltd.  All  rights  reserved, 
doi:  1 0.1 01 6/j.rser.2009. 07.01 1 


duration.  Sunshine  duration  can  be  easily  and  reliably  measured 
and  data  are  widely  available. 

The  design  of  a  solar  energy  conversion  system  needs  exact 
knowledge  regarding  the  availability  of  global  solar  radiation. 
Sunshine  hours  are  measured  at  many  locations  around  the 
world,  while  global  radiation  is  measured  at  selected  locations 
only.  In  order  to  overcome  this  defectiveness,  scientists  have 
developed  many  empirical  equations.  Most  of  the  sunshine 
based  this  equations  built  to  estimate  the  monthly  average  daily 
global  solar  radiation  are  of  the  modified  Angstrom-type 
equation. 

Determination  of  the  solar  energy  capacity  of  a  region  requires 
that  extensive  radiation  measurements  of  high  quality  be  made  at  a 
large  number  of  stations  covering  the  major  climatic  zones  of  the 
region.  In  this  regard,  recently,  several  empirical  formulas  using 
various  parameters  have  been  given  to  estimate  the  solar  radiation 
around  the  world.  Estimations  of  the  monthly  average  daily  global 
solar  radiation  for  a  large  number  of  locations  are  presented  in 
various  works  [1-63]. 

The  main  objective  of  this  study  is  to  review  the  global  solar 
radiation  models  available  in  the  literature,  including  the  study 
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Nomenclature 

a,  b  regression  coefficients 

H  monthly  average  daily  global  radiation  on  hor¬ 

izontal  surface  (MJ  itT2  d_1) 

H0  monthly  average  daily  extraterrestrial  radiation 
(MJ  m~2  d_1) 

Isc  solar  constant  (W  m~2) 

S  monthly  average  daily  hours  of  bright  sunshine  (h) 

S0  monthly  average  day  length  (h) 

Z  altitude  of  site  (km) 

<5  solar  declination  (°) 

9)  latitude  of  the  location  (°) 

a>s  sunset  hour  angle  (°) 


carried  out  on  the  estimation  of  the  monthly  average  daily  global 
solar  radiation  on  horizontal  surfaces. 


Model  1:  Glover  and  McCulloch  proposed  the  following 
equation  which  depends  on  9  and  is  valid  for  9)  <  60°  [4]: 

£-  =  0.29  cos(<p)  +  0.52  (jL)  (6) 

Model  2:  Page  has  given  the  coefficients  of  the  modified 
Angstrom-type  model,  which  is  believed  to  be  applicable  any¬ 
where  in  the  world,  as  the  following  [2]: 

k-M3+0-48©  |7> 

Model  3:  Rietveld  examined  several  published  values  of  the  a 
and  b  from  following  equations,  respectively  [5]: 

a  =  0.10  +  0.24^  (8a) 

b  =  0.38  +  0.08  (JPj  (8b) 

Model  4:  Iqbal  used  dada  obtained  from  three  locations  in 
Canada  to  propose  the  correlations  [6]: 


2.  Models  used 


The  simple  model  used  to  estimate  monthly  average  daily 
global  solar  radiation  on  horizontal  surface  is  the  modified 
form  of  the  Angstrom-type  equation.  The  original  Angstrom- 
type  regression  equation  related  monthly  average  daily  radia¬ 
tion  to  clear  day  radiation  at  the  location  in  question  and 
average  fraction  of  possible  sunshine  hours  [1],  Page  [2] 
and  others  have  modified  the  method  to  base  it  on  extra¬ 
terrestrial  radiation  on  horizontal  surface  rather  than  on  clear 
day  radiation  [64]: 


(1) 


where  H  is  the  monthly  average  daily  global  radiation,  Ha  is  the 
monthly  average  daily  extraterrestrial  radiation,  S  is  the  monthly 
average  daily  hours  of  bright  sunshine  (h),  SD  is  the  monthly 
average  day  length  (h),  and  a  and  b  are  empirical  coefficients. 

The  monthly  average  daily  extraterrestrial  radiation  on  a 
horizontal  surface  (H0)  can  be  computed  from  the  following  equation 
[64]: 


24 

Ho  =  —Ur  I  1  +0.033  cos 


360D\ 
~36lT  ) 


(  „  .  2jrtws  .  .  \ 

x  I  cos  9)  cos  a  sin  a>s  +  sin  <p  sin  6  I 


(2) 


where  Isc  is  the  solar  constant  (=1 353  W  m~2),  9)  the  latitude  of  the 
site,  <5  the  solar  declination,  cos  the  mean  sunrise  hour  angle  for  the 
given  month  and  D  the  number  of  days  of  the  year  starting  from 
first  January.  The  solar  declination  (5)  and  the  mean  sunrise  hour 
angle  (<os)  can  be  calculated  by  Eqs.  (3)  and  (4),  respectively  [64]. 


<5  = 


23.45  sin 


360(D  +  284)' 
365 


(3) 


u>s  =  cos  1  [-tan(<$)  tan(9>)[  (4) 

For  a  given  month,  the  maximum  possible  sunshine  duration 
(monthly  average  day  length,  S0 )  can  be  computed  by  using  the 
following  equation  [64]: 


So 


2 


15 


(5) 


The  regression  models  that  have  been  proposed  in  the  literature 
are  given  below. 


^  =  0.791  -0.635  (j^)  (9a) 

+  .0.163  +  0478(1) -0655g)2  (9b) 

Model  5:  Gariepy  has  reported  that  the  empirical  coefficient  a 
and  b  are  dependent  on  mean  air  temperature  (T)  and  the  amount 
of  precipitation  (P)  [7[: 

a  =  0.3791  -  0.0041  (T)  -  0.0176(P)  (10a) 

b  =  0.4810  +  0.0043(T)  +  0.0097(P)  (10b) 

Model  6:  Kilic  and  Ozturk  [8]  have  determined  that  the 
coefficients  a  and  b  are  a  function  of  the  solar  declination  (5)  in 
addition  to  both  9)  and  2,  as  given  by  the  equations  [11]: 

a  =  0.103  +  0.01 7Z  +  0.198  cos(<p-  8)  (11a) 

b  =  0.533  -  0.165 cos(ip-  8)  (lib) 

Model  7:  Lewis  obtained  the  following  linear  regression 
equation  to  estimate  the  daily  diffuse  radiation  for  three  stations 
in  Zimbabwe  [9[: 

^  =  0.754-0.654  (Jpj  (12) 


Model  8:  Kholagi  et  al.  derived  the  following  equations  from  the 
data  measured  at  three  different  stations  in  Yemen  [10]: 

jj-  =  0.191  +  0.571  (J^j  (13a) 


0.297 +  0.432  (JP\  (13b) 

Jt=  0.262  +  0.454^  (13c) 


Model  9:  Dogniaux  and  Lemoine  have  also  proposed  following 
equation,  where  the  regression  coefficients  a  and  b  seem  to  be  as  a 
function  of  cp  in  average  and  on  the  monthly  base  in  these 
equations,  respectively  [11]. 

a  =  0.37022 -0.00313(^)  (14a) 

b  =  0.32029  +  0.00506(9))  (14b) 
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H_ 

Ho 


=  (0.34507  -  0.00301$))  +  (0.34572 


0.00495$))  (jpj  for 


January 


H 

Ho 


H 

Ho 


H 

Ho 


H 

Ho 


=  (0.33459  -  0.00255$))  +  (0.35533 


+  0.00457 ip)  1  for  February 

=  (0.36690  -  0.00303$))  +  (0.36377 
+  0.00466$))  (Jpj  for  March 


=  (0.38557  -  0.00334$))  +  (0.35802 


0.00456 $))  J  for  April 


=  (0.35057  -  0.00245$>)  +  (0.33550 


0.00485$))  1  for  May 


=  (0.39890  -  0.00327$))  +  (0.27292 

HO 

+  0.00578$))  fj~)  for  June 


H 

Ho 


H 

Ho 


H 

Ho 


H 

Ho 


H 

Ho 


H 

Ho 


=  (0.41234  -  0.00369$))  +  (0.27004 
+  0.00568$))  (jpj  for  July 


=  (0.36243  -  0.00269$))  +  (0.33162 


0.00412$))  for  August 


=  (0.39467  -  0.00338$))  +  (0.27125 


0.00564$))  (  for  September 


=  (0.36213  -  0.00317$))  +  (0.31790 


0.00504$))  (  —  )  for  October 


=  (0.36680  -  0.00350$))  +  (0.31467 


0.00523$))  for  November 


=  (0.36262  -  0.00350$))  +  (0.30675 


0.00559$))  J  for  December 


(14c) 


(14d) 


(14e) 


(14f) 


(14g) 


(14h) 


(141) 


(14j) 


(141<) 


(141) 


(14m) 


(14n) 


Model  10:  Benson  et  al.  proposed  two  different  equations  into 
two  intervals  of  a  year  depending  on  the  climatic  parameters  [12]. 

=  0.18  +  O^of^  for  January-March  and 

1*0  \^0/ 


October-December 


(15a) 


^  =  0.24  +  0.53  for  April-September  (15b) 

i*o  \^oJ 


Model  11:  Ogelman  et  al.  have  correlated  (H/H0)  with  (S/SG)  in 
the  form  of  a  second  order  polynomial  equation  [13]: 


».  0,95  +  0.676(|) -0.1 42g 


(16) 


Model  12:  Ibrahim  obtained  the  following  equations  to  predict 
the  monthly  main  daily  diffuse  radiation  in  Cairo,  Egypt  [14]: 


£-0.79-0.59(1 


£  =  0.252  -  0.000,  (f)- 0.083  (|) 


(17a) 


(17b) 


Model  13:  Bahel  et  al.  suggested  the  following  relationship 
[15]: 


"-0.175  +  0.552(1 


(18) 


Model  14:  Zabara  proposed  monthly  a  and  b  values  of  the 
modified  Angstrom  model  as  a  third  order  function  of  maximum 
possible  sunshine  duration  (S)  and  day  length  (SG)  [16]. 

a  =  0.395  -  1.247^  +2.680 (jpj  -  1.674 (jpj  (19a) 


b  =  0.395  +  1.384 


:)  +2  055(s 


(19b) 


Model  15:  Jain  fitted  the  Angstrom  equation  using  the  least- 
square  method  to  the  monthly  average  daily  global  radiation  and 
the  sunshine  duration  data  of  31  Italian  locations.  The  equation 
using  the  mean  values  of  these  locations  is  given  as  follow  [17]: 


tj-  =  0.177  +  0.692 
”0  \^0 


(20) 


Model  16:  Alsaad  derived  the  Angstrom-type  equation  to 
estimate  the  monthly  average  daily  global  radiation  for  Amman, 
Jordan  [18]: 


«=0.,74  +  0.6,5(| 


(21) 


Model  17:  Bahel  also  developed  a  worldwide  correlation  based 
on  bright  sunshine  hours  and  global  radiation  data  of  48  stations 
around  the  world,  with  varied  meteorological  conditions  and  a 
wide  distribution  of  geographic  locations  [19]: 


"=0.,6  +  °.87(|) -06,(1)  +0.34(| 


(22) 


Model  18:  Gopinathan  has  suggested  that  the  coefficients  a  and 
b  are  a  function  of  (S/SD)  and  the  altitude  of  the  site  (Z),  as  given  by 
the  equations  following  equation  [20]: 


a  =  0.265  +  0.07Z-  0.135 


b  =  0.401  -0.108Z  + 0.325 


So 


(23a) 


(23b) 


Model  19:  Gopinathan  has  given  the  following  correlations 

[21]: 

=  -0.309  +  0.539  cos($>)  -  0.0693Z  +  0.290  (  f. 

llo  \Jo. 

1 .527  -  1 .027  cos ($>)  +  0.0926Z  -  0.359  (j)  j  (j 


(24) 


where  Z  is  altitude  in  kilometres. 
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Model  20:  Gopinathan  obtained  the  following  equations  from 
the  experimental  data  of  three  stations  in  India  [22]: 


^  =  0.931-0.814^ 
H  V^o. 


£-,.,94-0.83s(£)  -0.446(1 


(25a) 


(25b) 


Model  21 :  Jain  and  Jain  used  following  equation  to  estimate  the 
global  radiation  over  eight  Zambian  locations  [23]: 


£.0.240  +  0.513(1 


(26) 


Model  22:  Newland  suggested  following  equation  which 
includes  a  logarithmic  term  [24]. 


£.  0.34  +  0.40  (|)+ 0,7  log  (| 


(27) 


Model  23:  Soler  has  given  a  modified  Angstrom-type  equation 
for  each  month  as  follows  [25]: 


0  =  0.18  +  0.66^^- )  for  January 

Ho  \Jo 

0  =  0.20  +  0.60  )  for  February 

n0  \^oJ 

0-  =  0.22  +  0.58  )  for  March 

nn  V-Jo 


0  =  0.20  +  0.62  ( 0)  for  April 
no  \^oJ 

0  =  0.24  +  0.52  ( 0)  for  May 

no  \^n/ 


^  =  0.24  +  0.53  f|-')  for  June 
n0  \Jo 

~  =  0.23 +  0.53  0’')  for  July 
0  =  0.22  +  0.55  (0  )  for  August 

no  \^o 


0  =  0.20  +  0.59  f^-')  for  September 
no  \^oJ 

0  =  0.19 +  0.60  00  for  October 
no  \^oJ 

0  =  0.17  +  0.66(0)  for  November 
no  V-Jo 


0  =  0.18  +  0.65  rj’- )  for  December 

no  \^o 


(28a) 

(28b) 

(28c) 

(28d) 

(28e) 

(28f) 

(28g) 

(28h) 

(28i) 

(28J) 

(28k) 

(281) 


Also,  the  regression  coefficients  of  a  and  b  was  found  for  the  1 00 
stations  given  in  Ref.  [25]  as  follows: 


a  =  0.179  +  0.099 


b  =  0.1640  +  0.178600  -  1.09350 


(29m) 


(29n) 


Model  24:  Luhanga  and  Andringa  derived  their  own  model  as 
follow  [26]: 


0  —  0.241  +  0.488  rj*- 

n0  \50 


(30) 


Model  25:  Raja  and  Twidell  provided  the  following  equations 
using  the  data  from  five  main  observatories  in  Pakistan  and  by 
taking  into  account  the  effect  of  latitude,  tp,  [27,28], 


^  =  0.335  +  0.3670^ 


H 

Tfo 


=  0.388  cos(ip) 


°-367© 


(31a) 


(31b) 


Model  26:  Akinoglu  and  Ecevit  obtained  the  correlation  below 
between  (H/H0)  and  (n/N)  in  a  second  order  polynomial  equation 
for  Turkey  [29]: 


£-0.145  +  0.845(1) -0.280(1 


(32) 


Model  27:  Jain  explained  results  of  the  following  linear 
regression  analysis  of  measured  data  for  the  three  locations 
(Salisbury,  Bulawayo  and  Macerata,  Italy),  respectively  [30]: 


£.0.3,3 +  0.474g 
0  =  0.307  +  0.488  0^ 

0  —  0.309  +  0.599 
n0  \Jo 


(33a) 

(33b) 

(33c) 


Model  28:  Samuel  expressed  the  ratio  of  global  to  extra¬ 
terrestrial  radiation  as  a  function  of  the  ratio  of  sunshine  duration 
[31]: 


£.-0.14  +  2.52(1) -3.7, 


2.24 


(0 


(34) 


Model  29:  Tasdemiroglu  and  Sever  came  up  with  the  following 
equations  for  Turkey  in  general  [32]: 


0  =  0.622  -  0.35000  for  0.2  <  (JOj  <  0.94 


(35) 


Model  30:  Tasdemiroglu  and  Sever  also  developed  a  correlation 
between  (H/H0)  and  (S/SG)  in  a  second  order  polynomial  equation 
for  six  locations  (Ankara,  Antalya,  Diyarbakir,  Gebze,  Izmir  and 
Samsun)  of  Turkey  as  follow  [33]: 


£.0.225  + 0.0, 4g)  +0  00,(| 


(36) 


Model  31:  Louche  et  al.  suggested  the  model  below  to  estimate 
global  solar  radiation  [34]: 


^  =  0.206  +  0.546(0 

n0  \Jo 


(37) 


Model  32:  Lewis  derived  the  equations  including  the  linear  and 
three-order  polynomial  relationships  between  the  monthly 
average  values  of  (H/H0)  and  (S/SD)  for  locations  in  the  state  of 
Tennessee,  U.S.A.  as  follows  [35]: 


£  =  0.14  +  0.57(1 


£  =  0.8, -3.34(1) +7.3S(| 


4.51 


(38a) 


(38b) 


Model  33:  Gopinathan  and  Soler  suggested  linear  equation  for 
locations  with  latitudes  between  60  N  and  70°N  [36]: 


0  =  0.1538  +  0.7874(0 

n0  \-->o 

0  =  0.1961  +  0.72120 

n0  V50 


(39a) 


(39b) 


2584 


K.  Bakirci / Renewable  and  Sustainable  Energy  Reviews  13  (2009)  2580-2588 


Model  34:  Veeran  and  Kumar  obtained  the  following  linear 
relation  to  estimate  the  monthly  average  daily  global  radiation  at 
two  tropical  locations  (Madras  and  Kodaikanal,  India),  respectively 
[37]: 


tf- 0  34 +  0  32(1 


H 

Ho 


^  0.27  +  0.65(|) 


(40a) 

(40b) 


Model  35:  Yildiz  and  Oz,  using  the  measured  data  gathered 
from  five  stations  located  in  different  places  of  Turkey,  developed 
the  following  equation  [38]: 


=  0.2038  +  0.9236  (Jpj  -  0.391 1  ^ 


(41) 


Model  36:  Tiris  et  al.  also  suggested  the  following  correlations 
[39]: 


+  =  0.2262 +  0.41 8  + 
n0  V+o 


(42a) 


+  =  0.4177  + 0.0070  (+) -  1.9096  (+)  -1.19^)  (42b) 


^  =  0.1426-0.119^ 
+  =  0.0851  -0.298  (+ 

\Jo 

+  =  0.4428-0.1747  (+ 


(42c) 

(42d) 

(42e) 


Model  37:  Aksoy,  using  the  data  from  August  1993  to  July 
1995  obtained  from  the  Turkish  State  Meterological  Service, 
developed  a  quadratic  relationship  between  (H/H0)  and  (n/N) 
in  order  to  estimate  monthly  average  global  radiation  for 
Ankara,  Antalya,  Samsun,  Konya,  Urfa  and  Izmir,  Turkey,  as  follows 
[40]: 


A, 0.148  +  0.668  (i) -0.079(1) 


(43) 


Model  38:  Tiris  and  Tiris  derived  the  following  equations  from 
the  experimental  data  measured  in  Gebze,  Turkey  in  the  period 
from  January  1984  to  December  1992  [41]. 


^  =  0.652  -  0.482  ^  j  for  0.23  <  (jpj  <  0.76 


(44) 


Model  39:  Said  et  al.  obtained  the  following  equations  to 
estimate  monthly  average  global  and  diffuse  radiation  on  a 
horizontal  surface  at  Tripoli,  Libya  [42]: 

tj-  =  0.215  +  0.527 (+ 

Ho  Wo 


H 

Ho 

Hd 

H 


=  0.1+0.874 


=  1.625  -3.421 


Model  40:  Ampratwum  and  Dorvlo  suggested  the  following 
logarithmic  equation  for  Seeb  weather  station  in  Oman  [43]: 


+  =  0.6376  +  0.2490  log  ( + 
n0  V^o 


(46) 


Model  41 :  Togrul  and  Onat  developed  equations  to  estimate  the 
monthly  mean  global  solar  radiation  (H)  for  Elazig,  Turkey  [44]: 


H  =  -1.3876  +  0.518 H0  +  2.3064 


H  =  2.765  +  4.9597  sin(S)  +  2.2984 


(47a) 

(47b) 


Model  42:  Ulgen  and  Ozbalta  proposed  the  following  linear  and 
second  degree  equations  for  Izmir-Bornova,  Turkey  [45]: 


+  =  0.2424  +  0.5014  ( + 
H0  \^0 


U  =  0.0959  +  0.9958  (Jpj  -  0.3922  (J- 


(48a) 

(48b) 


Model  43:  Ertekin  and  Yaldiz  have  suggested  following 
polynomial  correlation  equations  for  Antalya  city  of  Turkey  [46]: 

w.  -  - 24 275  + 11 946 (I)  -  16-745(|)2  +  79575  (I)’ 

(49) 

Model  44:  Elagib  and  Mansell  have  investigated  the  possibility 
of  establishing  monthly-specific  equations  for  estimating  global 
solar  radiation  across  Sudan.  The  authors  have  reported  to  the  best 
performing  equations  for  each  month  as  given  following  [47]: 


+  =  0.1357  +  0.3204(<p)  +  0.0422(Z) 

‘‘O 


H 

Ho 


0.4947  ( — J  for  January 


=  0.1563  +  0.3166(<p)  +  0.1006(Z) 


0.4593  (!) 


for  February 


TJ  /  \  U.  /zbJ 

—  =  0.7727  (  —  for  March 

H o  \^oJ 


(50a) 

(50b) 

(50c) 


+  =  0.1640  +  0.0397(Z)  +  0.5773(+)  for  April  (50d) 

rl0  \^oJ 

+  =  0.0709  + 0.8967  (+)  -0.2258  for  May  (50e) 


H 

Ho 


-0.0348  +  1 .5078  Q  j  -  0.8246  j  for  June  (50f) 


+  =  0.3205  +  0.1444(^5)  +  0.0782 (h)  +  0.2916(" +  )  for  July 

n0  \^oJ 


H 


(45a) 

Ho 

\2 

0  255  fe 

) 

(45b) 

H 

Ho 

/ 

'  s\2 

+  2.1851 

.5) 

(45c) 

(50g) 


=  0.2720  +  0.0369(y?)  +  0.1 01 7(Z)  +  0.3888  (  for  August 


(50h) 


So 


=  -0.3710  +  2.57831  =- )  -  1.6788  (+=)  for  September 


(50i) 


^  =  0.1593  -  0.1043(99)  +  0.0609(Z)  +  0.5916(" +  )  for  October 
n0  \^oJ 

(50j) 
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£-  =  0.1786  +  0.0199(Z)  +  0.5441  (^-)  for  November  (50k) 

Hr,  \S„ 


H 

TTo 


=  0.1714  +  0. 1329((?)  +  0.0482(Z)  +  0.5015^) 


for  December 


(501) 


Model  45:  Chegaar  and  Chibani  have  suggested  two  models  for 
estimating  monthly  average  daily  global  on  a  horizontal  surface 
[48]: 


H 

Ho 

H 

Ho 

H 

TTo 


Model  46:  Ulgen  and  Hepbasli  developed  the  following  empi¬ 
rical  correlations  for  the  city  of  Izmir,  Turkey,  for  estimating  H  [49]. 


//  =  0.3092  cos  (<p)  +  0.4931  ( ^ 
Ho  wo. 


(52a) 


m  -  02854 + 02591  (!) + 0■6,7,  (!)  - 0  4834(!) 


(55b) 


Model  50:  Almorox  and  Hontoria  proposed  the  following 
exponential  equation  [53]: 


^-=  -0.0271  +0.3096  exp  ( ^ 
Ho  Vyo 


(56) 


Model  51 :  Ahmad  and  Ulfat  have  suggested  to  first  and  second 
order  polynomial  equations  developed  for  Karachi  of  Pakistan. 


for  Algiers  and  Oran 

(51a) 

The 

H 

for  Beni  Abbas 

(51b) 

Ho 

H 

for  Tamanrasset 

(51c) 

Ho 

-0.405 


=  0.348  +  0.320 


(57a) 

0.070  (pj  2  (57b) 


Model  52:  Rensheng  et  al.  suggested  a  logarithmic  relationship 
between  the  daily  global  radiation/daily  extraterrestrial  solar 
radiation  (H/H0)  and  the  temperature  difference  between  the 
maximum  and  minimum  daily  air  temperature  ( TM  -  Tm )  as 
equations  given  following  [55]: 


£  -  0.2408  +  0.3625  (|)  +  0.4507  (|) '  -  0.3708  (f)  ’ 

(52b) 

Model  47:  Togrul  and  Togrul  suggested  the  following  equations 
which  they  obtained  from  the  relation  between  S/SG  and  H/H0  by 
trying  different  regression  types  for  Ankara,  Antalya,  Izmir, 
Yenihisar  (Aydin),  Yumurtahk  (Adana)  and  Elazig  in  Turkey  [50]: 

^  =  0.318  + 0.449  (53a) 

+^  =  0.698 +  0.2022  In f(T)  (53b) 

Ho  \^o J 

«-  0.1541  +  1.1714(|)  -0.705(|) 2  (53c) 

W.  ~  0.1796  +  0.98,3(1)  -  0.2958(|)!  -  OOSSog)’ 

(53d) 

//  =  0.3396  e0  S985(S/S„)  (53e) 

Ho 

it  /  r  \  0.4146 

Hr°-7316u)  (53f) 

Model  48:  Akpabio  and  Etuk  suggested  an  Angstrom-type 

correlation  equation  given  as  following  by  using  measurements  of 
global  solar  radiation  and  sunshine  duration  data  during  the  period 
from  1984  to  1999  at  Onne  (within  the  rainforest  climatic  zone  of 
southern  Nigeria)  [51]: 

£.0.23  +  0  38(1)  (54) 

Model  49:  Ulgen  and  Hepbasli  suggested  the  following  the 
linear  and  third  order  polynomial  equations  for  Ankara,  Istanbul 
and  Izmir  in  Turkey  [52]: 

-^-  =  0.2671  +0.4754^  (55a) 

Hq  \^oJ 


^-  =  a\n(TM-Tm)  +  b(£jC +  d  (58) 

where  a,  b,  c  and  d  were  empirical  coefficients  which  they  gave 
their  study  [61], 

Model  53:  Almorox  et  al.  reported  the  monthly-specific 
equations  for  estimating  global  solar  radiation  from  sunshine 
hours  for  Toledo,  Spain  as  given  following  [56]: 


=  0.285  + 0.444 for  January,  (59a) 

Ho  \^oJ 

„  =  0.272  +  0.465 ( §-]  for  February  (59b) 

Ho  \$0  J 

=  0.291  +  0.491  (-£)  for  March  (59c) 

Ho  \^oJ 

=  0.266  +  0.495  ( ^  )  for  April  (59d) 

H0  \^o  / 

£  =  0.286  +  0.475 ( §-]  for  May  (59e) 

Ho  \^0 J 

^-=0.311+0.439^)  for  June  (59f) 

Ho  \^oJ 

=  0.329  +  0.406  (jpj  for  July  (59g) 

^1=0.313  +  0.410^)  for  August  (59h) 

Ho  wo / 

^  =  0.271  +0.479 f for  September  (59i) 

Ho  \$0  J 

=  0.259  +  0.465  ( -+)  for  October  (59j) 

Ho  \^oJ 

yj-  =  0.279  +  0.431  (^7+)  for  November  (59k) 

H0  \^oJ 

5  =0.282  +  0.428  (£-)  for  December  (591) 

H0  \-->0  J 
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Model  54:  Tahran  and  Sari  have  suggested  two  models  to  predict 
solar  radiation  over  the  Central  Black  Sea  Region  of  Turkey.  These 
quadratic  and  cubic  polynomial  models  are  given  as  following, 
respectively  [57]: 


^-  =  0.1874  +  0.8592  (jpj  -0.4764^- 


(60a) 


©■11126©i+045l6©>  |60b> 


^  =  0.1520  +  1.13341  ^  )  -  1.11261  ^ 

Model  55:  Jin  et  al.  have  suggested,  based  on  the  radiation  data 
and  the  geographical  information  including  latitude  and  altitude  at 
all  69  stations  in  China,  nine  general  equations  are  given  as 
following  [58]: 


^  =  0.1332  + 0.6471 
Ho  Vpo 


^  =  0.1404  +  0.6126^ 

Ho  \3o 


0.0351 


^  =  0.1275  +  0.7251  ( )  -  0.2299 
Ho  \50 


H 


©’ 


0.1837 


So 


„  =  0.0855  +  0.0020(0))  +  0.030Z  +  0.5654. 

Ho  \yo 


^  =  2.1186  -  2.0014 cos(yj)  +  0.0304Z+  0.5622 ( $ 
Ho  V^O. 


(61a) 

(61b) 
(61c) 
(61  d) 
(Gle) 


Model  57:  Rensheng  et  al.  have  suggested  new  equations, 
basing  on  the  Angstrom  model  and  the  Bahel  model,  by  using  the 
daily  global  radiation  data  and  sunshine  hours  from  1994  to 
1 998  at  86  stations  in  China.  These  equations  are  given  as  following 
[60]: 


ip-  =  0.176  +  0.563 

H0  V-JO 


(63a) 


=  (0.122  +  0.001(0))  +  2.57  x  10~2Z)  +  0.543  ( f]  (63b) 

Ho  \^o/ 


H 

Ho 

H 

Ho 


=  (0.280  -  0.141  cos (<p)  +  2.60  x  10~2Z)  +  0.542  Q  'j  (63c) 


=  (0.275  +  4.27  x  10~5X  -  0.141  cos (<p)  +  2.63 


x  10~2Z)  +  0.542^ 

WO 


^  =  (0.117  +  4.11  x  10~57.  — O.OOl0)+2.59  x  10~2Z)  +  0.543 .  _ 
H0  \yo 


(63d) 

S\ 

So 

(63e) 


^-=  [-0.196  cos{(p)  +  2.2  x  10~2Z  +  0.329] 
Hq 


[0.097  cos(0>)  +  6.72  x  10~3Z  + 


04571  © 


(63f) 


H 

Tfn 


=  (0.1094  +  0.0014(01)  +  0.0212Z)  +  (0.5176  +  0.0012(0)) 


0.01 50Z) 


© 


(61  f) 


=  (1.8790  -  1.7516cos(0j)  +  0.0205Z)  +  (1.0819 
Hn 


-  0.5409  cos(<p)  +  0.0169Z)  I 


(61  g) 


=  (0.0218  +  0.0033(0?)  +  0.0443Z)  +  (0.9979  -  0.0092(0?) 
Ho 


>© 


-  0.0852Z)  —  +  (-0.5579  +  0.0120 (<p)  +  0.1005Z)  — 


© 


(61  h) 


H 

Ho 


=  (4.2510  —  4.1878cos(0))  +  0.0437Z) 

+  (-10.5774  +  1 1.4512  cos(0>)  -  0.0832Z) 


(12.7247  -  13.0994  cos(0?)  +  0.1000Z) 


So 


(61  i) 


where  H  is  the  monthly  average  daily  global  radiation,  Ha  is  the 
monthly  average  daily  extraterrestrial  radiation,  S  is  the  actual 
sunshine  duration,  S„  is  the  maximum  possible  sunshine  duration, 
tp  is  the  latitude  of  the  site,  Z  is  the  altitude  of  the  site  in  kilometers. 

Model  56:  Aras  et  al.  have  suggested  following  linear  and 
polynomial  correlation  equations  to  use  generally  for  twelve 
provinces  in  the  Central  Anatolia  Region  of  Turkey  [59]: 


^  =  0.3078  +  0.4166^ 
H0  \Jo 

-tj-  =  0.3398  +  0.2868  ( ^ 
Ho  \^o 

^-  =  0.4832-0.6161^ 


0.1187 


(62a) 

(62b) 


1 8932 ©2 0975 ©’  |62c> 


H 

TT0 


-  [0.000U- 0.195 cos(0>)  +  2.28  x  IO^Z  +  O.313] 
[—0.00027.  +  0.097  cos(0))  +  5.69  x  10~3Z  +  0.476]  (Jpj 


H 

H, 


—  =  [0.000U  +  0.0020)  +  2.27  x  10~2Z  +  0.094] 


+  [—0.00027.  —  0.00080)+  5.36  x  10  3Z  +  0.586] 

\3o 


H 

H„ 

H 

TTo 


H 

TT0 


H 

TTo 


^,50+,i45©-,  474(£)  +o-963(s 

-  (0.0010)  +  2.41  X  10~2Z  +  0.109)  +  1.029^ 

-,.2,6©2+0.7S7(| 


[-O.112cos(0)) +  2.43  x  10~2Z  +  0.234]  +  1.026 


-1.2°9(|)  +0.782(| 


=  (0.0007 0)  +  2.44  x  10~2Z  -  0.0057.  +  2.24  x  lO^5© 


■0.370)  +  ,. 026(1) -1208  ©2  + 0.783(1)’ 


(63g) 


(63h) 

(631) 

(63j) 

(63k) 

(631) 


=  [-0.087  cos(0))  +  2.44  x  10~2Z  -  0.004A  +  1.86  x  10~5Z2 
Ho 

+  0.426]  +  1 .024  (Jpj  -  1 .204  (Jpj  +  0.779  ^  (63m) 


K.  Bakirci / Renewable  and  Sustainable  Energy  Reviews  13  (2009)  2580-2588 


2587 


H_ 

TTo 


[—0.233  cos((p)  +  2.64  x  1CT2Z  +  0.336] 

+  [2.140  cos (<p)  -0.1  Z  -  0.670]  ^  +  [-5  cos (<p)  +  0.3 Z 
-  2.744]  (-S-)  +  [3.042  cos(tp)  —  0.2Z  —  1 .638]  (S-\  (63n) 


K-  =  [~0.141cos(p)  +  2.84  x  1(T2Z  -  0.014A.  +  6.7  x  1(T5Z2 
H0 

+  1 .01 2]  +  [1 .740  cos ((P)  -0.1  Z  +  0.0697.  -  0.003X2 

-  4.061]  (J- J  +  [-3.867  cos(<p)  +  0.3Z  -  0.1997.  +  0.0009X2 

+  12.402]  (^-)  +  [2.115cos(^)  -  0.2Z  +  0.1647.  -  0.0008Z2 

-  9.442]  (0 3  (63o) 

where  X,  tp  and  Z  (in  kilometers)  are  the  longitude,  latitude  and 
altitude  of  the  used  stations. 

Model  58:  Bakirci  has  suggested  following  polynomial  correla¬ 
tion  equations  for  Erzurum  city  of  Turkey  [61]: 

^  =  0.6307  -0.7251^0 +1.2089^  -0.4633(0)  (64) 

Model  59:  Bakirci  reported  equations  for  estimating  global  solar 
radiation  from  sunshine  hours  for  Erzurum,  Turkey  as  follows  [62]: 

=  0.6716  +  0.0760nT)  (65a) 

Hq  V-JQ  / 


^  =  0.5622  +  0.5444  0j  -  0.4490  0j  (65b) 

Model  60:  Bakirci  reported  the  original  Angstrom-type 
equations  including  the  linear,  second-order  and  fifth-order 
polynomial  relationships  between  the  monthly  average  values 
of  (H/Hc)  and  (S/SG)  as  follows  [63]: 

-^  =  0.7836  -  0.0460  ( )  (66a) 

He  \^oJ 


~  =  1 .0192  -  1 .0547  0j  +  0.9661  0  j  (66b) 

«,-,1.225  +  ,  280,0(|) -5,  S.900(|)2 

+  994.7300)  -  920.350  0)  +329.93  0)  (66c) 

3.  Results  and  discussion 

In  this  study,  the  regression  constants  have  been  generally 
computed  using  observations  of  sunshine  hours  and  monthly 
average  daily  global  radiation  in  given  a  location.  In  the  studies 
carried  out  in  this  field  by  many  authors,  the  models  used  to  estimate 
the  monthly  average  daily  global  solar  radiation  on  a  horizontal 
surface  can  be  categorized  in  four  groups,  as  seen  from  Eqs.  (6)  to 
(66).  These  groups  consist  of  the  relations  derived  from  the 
Angstrom-type  equations  and,  the  types  of  the  models  are  given 
as  follows: 

Group  I  (linear  models):  The  models  derived  from  the  Angstrom- 
type  regression  equation  was  called  the  linear  models  because 


the  empirical  coefficients  a  and  b  were  obtained  from  the 
results  of  the  first  order  regression  analysis  (such  as  Model  2, 8, 
13  and  16).  (ii)  Group  II  (polynomial  models):  Some  researcher 
suggested  that  the  modified  Angstrom  type  relation  is  a 
second,  three  and  bigger  order  polynomial  equation  to 
estimate  the  monthly  average  daily  global  radiation  on  a 
horizontal  surface  (such  as  Model  28, 30,  54  and  60).  (iii)  Group 
III  (angular  models):  There  are  the  angular  models  derived  by 
modifying  the  original  Angstrom-type  equation  (such  as  Model 
1 , 6, 9  and  19).  (iv)  Group  IV  (other  models):  Special  cases  of  the 
modified  Angstrom-type  equation  were  categorized  in  this 
group.  These  cases  include  a  logarithmic  term,  non-linear 
model  and  exponential  equation  (such  as  Model  47,  52,  55 
and  57). 

4.  Conclusion 

Solar  radiation  data  are  essential  in  the  design  and  study  of  solar 
energy  conservation  devices.  In  this  regard,  empirical  correlations 
are  developed  to  estimate  the  monthly  average  daily  global 
radiation  on  a  horizontal  surface.  Sunshine  based  models  are 
employed  for  estimation  global  solar  radiation  for  a  location.  The 
correlation  equations  given  in  this  study  will  enable  the  solar 
energy  researcher  to  use  the  estimated  data  with  trust  because  of 
its  fine  agreement  with  the  observed  data. 

The  following  conclusions  may  be  drawn  from  the  present  study: 

(1)  Solar  energy  technologies  offer  a  clean,  renewable  and 
domestic  energy  source,  and  are  essential  components  of  a 
sustainable  energy  future.  In  the  design  and  evaluation  of  solar 
energy,  information  on  solar  radiation  and  its  components  at  a 
given  location  is  needed.  In  this  regard,  solar  radiation  models 
are  of  big  importance. 

(2)  Most  of  solar  radiation  models  given  to  estimate  the  monthly 
average  daily  global  solar  radiation  are  of  the  modified 
Angstrom-type  equation. 

(3)  It  may  be  concluded  that  the  models  presented  in  this  study 
may  be  used  reasonably  well  for  estimating  the  solar  radiation 
at  a  given  location  and  possibly  in  elsewhere  with  similar 
climatic  conditions. 
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